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(54) A digital time-delay beamformer 
for sonar systems 

(57) In a digital sonar receiver having a 
beamformer 38, own doppler nullifier 
(O.D.N.) 48, and adaptive time varied 
gain controller 36, the beamformer in- 
terrogates a multi-element transducer 
by utilizing an interleaved sampling 
scheme which samples the elements in 
the transducer array sequentially, in 
such a fashion that continuous, uni- 
formly spaced samples are obtained, 
converted into digital format, weighed 
and summed. The O.D.N. 48 offsets the 
frequency of the beam output signal 



and a digital generated signal represen- 
tative of the own ship's doppler and 
multiplies these two offsetted signals 
resulting in a difference frequency sig- 
nal corresponding to the target doppler. 
The adaptive time varied gain control 
apparatus 36 initially sets the gain of 
the preamplifiers 32 to be inversely 
proportional to a standard reverbera- 
tion curve during the first ping-cycle 
and subsequently updates the gain of 
the preamplifiers during each succes- 
sive ping-cycle by computing a gain 
correction from the average value of a 
predetermined number of beam output 
signals. 
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SPECIFICATION 



A digital time-delay beamformer for sonar systems 

5 This invention relates to a digital sonar receiver comprising a digital beamformer, a digital own doppler 5 
nullifier (O.D.N. ). Th sonar receiver receives sonar signals in analogue form and converts the sonar signals 
into digital form in the beamformer when forming beam output signals. The remainder of the sonar receiver 
which includes the O.D.N, and adaptive time varied gain controller operates digitally. 
In the past, sonar receivers for the most part have processed sonar signals in analogue format. However, 
10 more recent sonar receivers have been directed towards processing the signals in digital format. 10 
Previous digital beamformer designs have been directed towards a method of forming simultaneous 
beams. To form simultaneous beams a cylindrical array of staves had been simultaneously sampled at some 
initial time. The staves are then sampled simultaneously at later moments in time corresponding to the delay 
times introduced by the geometry of the array. This method of simultaneous beamforming when 
15 implemented, requires for each stave a sample and hold circuit and an AID converter. This method of 15 
beamforming has resulted in rather uneconomical use of time as well as hardware. 

The received sonar signals prior to beamforming comprise ambient noise and noise from any undesirable 
targets. The received sonar signal besides including an echo from the desired target also includes noise due 
to reverberation. Reverberation is a function of the range or time taken for a transmitted signal to hit a target 
20 and return. Reverberation decreases as range or time increases. Eventually the noise due to reverberation 20 
decreases below the ambient noise in the medium of operation which then becomes the dominant noise. 

To eliminate the effect of noise due to reverberation the gain of the preamplifiers in a sonar receiver is 
adjusted to be inversely proportional to an assumed "standard" reververation curve. The adjustment of the 
gain of the preamplifiers with respect to the "standard" reverberation curve has not been too effective in 
25 reducing the noise due to the actual reververation. Large discrepancies between the "standard" 25 
reverberation curve and actual reverberation easily result from surface reflections, back scattering, and 
marine life. 

Subsequent to beamforming, beam output signals that contain information corresponding to the target 
are often shifted in frequency. This frequency shift is known as the total doppler shift which, for the most 

30 part, is due to the doppler shift produced by target movement and the doppler shift produced by the source 30 
of transmission's movement or the own ship's doppler. It is necessary to eliminate the doppler effect of the 
own ship from the beam output signal so that the target information may be analyzed. Analogue sonar 
receivers eliminate the own ship's doppler by mixing the beam output signal with a generated signal 
corresponding to the own ship's doppler. The signal corresponding the own ship's doppler may be 

35 calculated from the velocity of sound in water, the ship's speed and heading which are all known. The result 35 
of this mixing process yields a sum signal and a difference signal where the difference signal is 
representative of the target information only. However, when the beam output signals are in digital form the 
effects of noise folding over from the negative frequency band may introduce error in the digital difference 
signal. Moreover the sum digital signal may overlap or lie too close to the difference digital signal thereby 

40 preventing the economical use of filters to filter out the sum signal. 40 
The present invention relates to a digital sonar receiver having a novel beamformer for digitally 
beamforming output signals from received sonar signals, a novel own doppler nullifier, and an apparatus for 
varying the gain of sonar preamplifiers to be inversely proportonal to the momentary average reverberation 
level of the medium in which the sonar receiver operates. 

45 in accordance with one aspect of this invention there is provided a beamformer for forming a plurality of 45 
digital beam output signals in a sonar receiver which has a multi-element symmetrical transducer array. 
Each beam output signal has a direction which is related to at least one respective directional element in the 
transducer array. The total number of samples to be taken is given as the product of the total number of 
beams M to be formed and the number sample-sets n desired to form any beam. Due to the symmetry of the 

50 array a sample-set may consist of two samples from each element symmetrically disposed on opposite sides 50 
of the at least one directional element. The beamformer sequentially samples the elements, in such a fashion 
that continuous, uniformly spaced samples are obtained. To obtain uniformly spaced samples the method of 
beamforming and beamformer apparatus sample each i element symmetrically disposed to the left and right 
of the directional element for each beam m during each consecutive sampling interval, z, where i is defined 

55 as an integer value greater than or equal to 0 and less than or equal to n-1 . The sampling interval z ts 
defined asan integer O^z^Z. The beamformer determines from the sampling interval which elements 
symmetrically disposed to the left and right of the directional element are to be sampled. The directional 
element which is utilized to form the beam m may be given by: 



55 



60 

m=Remainder 



60 



+M 

n 

M 



+ 1 




2 GB 2 029 574 A 



where di is a time delay of the i th sample to the right and left of beam m with respect to its first sample. For 
any given sampling interval beam m may be determined if the time delay is a funct.on of the sampling 
interval. The time delay is an approximation of the exact time delay for the movement of a sonar s.gnal 
through the medium of operation to imping on the elements of the sample-set after the sonar s.gnal has 
5 impinged on the at least one respective directional element for the beam. It should be understood that for 
this method the time delay is expressed as an integer number of sampling intervals where d. is chosen to be 
the integer number of sampling intervals closest to the exact time delay such that, the 
when divided by the total number of sample-sets to betaken has a different integer value fo each value of i. 
The value of i may be determined from the remainder [di/n] which is equal to the remainder [z/nl. During 
10 each sampling interval the beamformer sums together the samples taken from each i element to the right 
and left of the beam m. so as to form a sample-set. The beamformer converts the sample-set from analogue 
format to digital format. The beamformer multiples the digital sample-set with weights The b *»™™™* r 
adds each digital sample-set to its corresponding partial sum of sample-sets and stores this partial sum. The 
output of the partial sum for each beam when all n samples of the sample-set for beam m have been 
15 summed is the beam output signal corresponding to beam m. f «ui„„, h!1 . mnt i„„ rai , 
It should be understood that the sampling of the beams sequentially, m such a fashion that continuous, 
uniformly spaced samples are obtained may be determined for any multi-elements symmetrical array and 
mav be stored as a timing schedule in a memory device. It should therefore become apparent that a 
bTrnformer using !hTs sampling method only requires one A/D converter, .t should also be apparent that the 
20 sampling method results in an economical use of equipment and time. _ 

In accordance with another aspect of this invention there is provided an adapt.ve time varying gam control 
apparatus for varying the gain of the preamplifiers of the sonar receiver during each consecutive pmg-cycle 
to be inversely proportional to the momentary average reverberation in the medium of operation. The 
apparatus comprises a storage means which has a plurality of storage spaces for storing a stored gam. A 
25 preset means is provided for initially storing in these storage spaces gam step approximations of an inverse 
standard reverberation curve where each of the gain step approximates has a respective t.me duration which 
is determined from the standard reverberation curve. Also included in the apparatus is a means for storing 
each time duration and for associating each of the time durations with the storage space where the t me 
durations respective gain step approximation has been initially stored. A control means is provided for 
30 varying the ^gaTn of thf preampliSers to correspond to each stored gain during each ping-cycle wherebyeach 
stored gain is applied to the preamplifiers for a duration which is equal to the t.me duration associated with 
the stored gains respective storage space. An averaging means is also provided in the apparatus for 
averaging a plurality of beam output signals from a beamformer in the sonar rece.ver dunng each t.me 
duration The apparatus includes a means for determining a gain correction as a function of " reference 
35 signal and the average of the beam output signals during each time duration The apparatus further includes 
means for computing a new gain by summing the gain correction dur.ng each t.me duration with the 
respective stored gain associated with that time duration. Lastly, a means for storing a new gam in the stored 
gains respective storage space during each time duration is provided. Alc „ tho 
It should be understood that the gain step approximations may be umform step approximations. Also the 
40 means for determining a gain correction may be a means for computing a logarithmic funct.on of the 
Average of the beam output signal when divided by the reference signal. The determining means may also 
include a comparing means which compares the average of the beam output signals with reference s.gnal to 
give a coded value of the gain correction. This coded value may be used to access a memory which has 
precalculated values for the gain correction stored therein. .... .,-^ 0 w,, p:ri „ 45 

45 The adaptive time varying gain control apparatus continually updates the gam of the preamplifiers during 
each ping-cycle by correcting during each time duraction the gain approximations of the standard 
reverberation curve on the subsequently stored gains when the gain approx.mat.ons or the subsequently 
stored gains are applied to the preamplifiers. This adaptive time varied control apparatus eliminates no.se 
due to the momentary average reverberation in the medium of operation. llFarfnr 50 

50 In accordance with another aspect of this invention, there is provided an own doppler nullifier for 

extracting the own ship's doppler information from a digital beam output s.gnal of a beamformer to obtain 
target doppler information. The own doppler nullifier introduces an offset frequency A to the d g.tal beam 
output signal which offset is a function of the beamformer sampling frequency. The own doppler null.f.er 
apparatus and method also provides for the generation of a digital signal representative of the own ship s 
55 doppler information and offsetting the digital signal by a frequency amount B which is a funct.on of the 
beamformer sampling frequency and offset frequency A. The generated digital signal and the beam output 
signal are multiplied together to obtain a signal having a difference frequency component and a sum 
frequency component. Both frequency components will lie in their respective d.ff rence frequency and sum 
frequency bands. The difference frequency band will have a center frequency which is offset by an amount 
60 A-B whereby the affect of noise due to fold-over from the negative frequency band into the difference 
frequency band is eliminated. The own doppler nullifier provides for the sum frequency component to be 
filtered out leaving only the difference frequency component corresponding to the target doppler 

'"ThTowndoppler nullifier makes it possible to digitally filter the sum frequency component and difference 
65 frequency component of the mixed generated digital signal and digital beam output s.gnal. The introduction 
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of the offsets shifts the difference frequency component away from any affects of noise due to foid-over from 
the negative frequency band. Moreover, the offsets introduced to the signals shifts the sum frequency 
component further away from the difference frequency component to make the filtering of the sum 
frequency component from the difference frequency component possible. 
5 It should be understood that the offset frequency A may be chosen as the difference between the 5 
transmitted frequency of the sonar receiver minus the nearest lower multiple of the beamformer sampling 
frequency to that transmitted frequency. It should be understood that offset frequency B may be chosen to 
be equal to half the beamformer sampling frequency less offset frequency A so as to provide the optimum 
result. 

10 For a better understanding of the nature and objects of the invention, reference may be had to the 10 
following detailed description taken in conjunction with the accompanying diagrammatic drawings in 
which: 

Figure 1 shows a block diagram for a digital sonar receiver; 

Figure 2 is a schematic representation of a circular transducer array showing a wave front impinging the 
15 transducer array; 15 
Figure 3 is a timing diagram showing a partial sampling schedule for the digital beamformer; 
Figure 4 shows a schematic diagram of the digital beamformer; 

Figure 5 shows a "standard" reverberation curve, and the affect of ambient noise on such a curve (on the 
same sheet as Figure 2); 

20 Figure 6 is a block diagram for the input section of the sonar receiver; 20 
Figure 7 is a block diagram for the adaptive time varied gain controller; 

Figures 8A and 8B show the affect of noise on the difference band before and after the band is shifted by 
frequency amount A-B; 
Figure 9 is a block diagram for the own doppler nullification apparatus. 

25 Referring now to Figure 1 , a block diagram for the digital sonar receiver 10 is shown. The digital sonar 25 
rece i ver 70 receives sonar signals from a multi-element transducer array corresponding stave line 1 2. The 
stave line 12 apply the sonar signals to the input section 20 of the receiver 10. Input section 20 reduces the 
dynamic range of the signals, converts the signals to digital and forms beam outputs for subsequent 
analysis. Processing section 25 analyses the digital beam output data and extracts target information is then 

30 fed to output section 30 and subsequently onto a display or any other suitable means, for example, a 30 
recorder. 

The input section 20 of receiver 10 comprises a plurality of staves in a circular multi-element transducer 
array (not shown) where each stave is connected to a preamplifier in preamplifier circuit 32 by line 1 2. Each 
output of the preamplifier is connected to the input of an anti-aliasing filter of filter circuit 34. In order to 

35 reduce the large dynamic range at the input of receiver 10, each preamplifier of circuit 32 has a variable gain 35 
which is, by means of an adaptive time-varied gain controller 36, adjusted to be inversely proportional to the 
momentary value of the average reverberation level. The adjustment of the preamplifier circuit's gain 
achieves reverberation background level normalization. The time varied gain controller receives information 
via the beamformer 38, which samples the received signals, as to the momentary reverberation level. 

40 Reducing the dynamic range of the sonar signals reduces the digital word length required to represent the 40 
signal in the subsequent processing. 

The sonar signals once having their background noise level normalized by preamplifier circuit 32 are fed to 
the anti-aliasing filter circuit 34. Circuit 34 comprises a plurality of anti-aliasing filters where for each 
preamplifier there is a corresponding filter. Each anti-aliasing filter includes a band pass filter, with the center 

45 frequency set to correspond to the frequency of the transmitted signals. The width of the band frequencies 45 
passed by the filter will be determined by the maximum expected doppler shift between the target echo and 
the transmitted signal. The function of the anti-alias filters is to eliminate the passing of signals in the so 
called "alias" frequency bands, introduced by the subsequent sampling process. At the same time the signal 
to noise ratio is improved, because of the reduction in the bandwidth. As a result of the reduction in the 

50 dynamic range the sampling and subsequent digital processing can be performed more economically. 50 
Beamformer 38 samples the signals through multiplexers 40 and converts the signal into digital format by 
A/D converter and summer 42. The beamformer forms one beam for each related stave by sampling the 
staves symmetrically disposed about the related stave as well as the related stave itself and weighting the 
samples before adding them together to form the beam. Timing and control circuit 44 controls sampling by 

55 implementing a sampling method, which interrogates the staves in an interleaved fashion and forms the 55 
beams sequentially. 

The full beam output signal, in digital format, is passed through automatic gain controller 46 of processing 
section 25. The automatic gain controller increases or decreases the gain of each beam independent of the 
other beams after a comparison of the actual beam output with the average output of a number of preceding 
60 samples for that beam. 6 ^ 
The gain adjusted beam output, is passed through own doppler nullifier device 48 which eliminates the 
effect of the own ship's doppler on the frequency of the signal. 

The signal now enters a signal processor 50 which correlates the received signal with a plurality of 
predetermined independently stored signals, so called "replicas". The best correlation will be indicative of 
65 the velocity of the target in that beam. Such a correlation is performed for each beam in sequence. The 65 
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outDUt is sent to the output processor 52 of output section 30. 

Afine bear £ compufer 54 is included in processor section 25, which combines information from the 
^miS.TSStiS 46 and the signal processor 50 to more accurately compute « . direction 
within a beam the target information is received. This information ,s fed into the output : p ^*?°' 52 - 
5 The outout processor 52 receives information regarding the beam, target velocity and time since 
transmission and combines this information with inputs indicating the velocity of sound , m a "d the 

shTp's course From this the output processor 52 produces a digital signal which « 'nd.catrve of the range. 

upon the transducer. This delay may be expressed as: 

d=rn-cos (e-a))10 +6 ^iS, 
c 

where: r = the radius of the transducer in feet; 

i — woioritv of sound in water in feet/second; 
20 6 = the aZebZeen the direction of travel of the incident wavefront and the axis of the a rrav On 

a = Wangle between the axis of the stave considered and the axis of the array (in degrees). 

25 w^TJ^^oi axis of the array ,6.=0°> will generate signals in the 13 staves, delayed as 
shown in Table 1 . 



30 

Stave No. 

S7 

S6 
35 S5 

S4 

S3 

S2 

S1 
40 S36 

S35 

S34 

S33 

S32 
45 S31 



Location 
in array 

C+6 
C+5 
C+4 
C+3 
C+2 
C+1 

C(center 

C-1 

C-2 

C-3 

C-4 

C-5 

C-6 



Tabie I 



a (degrees) 

+60 
+50 
+40 
+30 
+20 
+ 10 
0 

-10 
-20 
-30 
-40 
-50 
-60 



Delay in uS. 

204.5 
146.0 
95.7 
54.8 
24.7 
6.2 
0 
6.2 
24.7 
54.8 
95.7 
146.0 
204.5 
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Bv takinq samples of the individual stave signals at moments in time corresponding to the delay times for 

'^hTbeamformer 38 of Figure 1 utilizes a method of interleaved sampling. The interleaved sampling 
mlthod sampleT^divrduafstave elements at different moments in time which approximate the de lay . mes 
55 of Table !.The stave elements are sequentially sampled, in such a fashion that continuous, uniformly spaced 

To formre^mcorresponding to each stave element of Figure 2, 7 samp.e-sets must be taken for each 
beam. Therefore the total number of samples to be taken is: 



60 



Z = Mn 
= <36)<7) 
= 252 



where M is th total number of beams to be formed (equal to 36) and n is the total number of san 
65 EJS to 7) Because each sample is taken at one moment in time, Z represents the total numb, „, moments 
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in time samples are to be taken. 

The sampling method samples each i element symmetrically disposed to the right and left of the beam m, 
being formed during each consecutive sampling interval z. It should be understood that o i n — 1 (or o i 6 for 
n = 7). Each sampling interval corresponds to one moment in time of the 252 different moments in time for 
5Z. 

The time duration of each sampling interval is chosen such that the received sonar signal may be 
adequately reproduced digitally. A total bandwidth deviation of 800 Hz for the received signal from the 
transmitted signal requires a sampling frequency of at least 1 ,600 Hz to adequately reproduce the received 
signal digitally. The time duration of the sampling interval is chosen to be a function of the first delay of 6.2 
10 f^S from Table I. Table II shows the calculation of different time durations. The inverse of the time duration is 
also given which represents the sampling frequency. From Table II a sampling frequency of fs=3200 Hz is 
chosen giving the time duration of each sampling interval as 1 .24 uS. 



10 



Table II 



15 



20 



25 



Time Duration of the 
Sampling Interval 
xo= 6.2 iiS 
X 

6.2 

3.1 

2.07 

1.55 

1.24 

1.04 

0.89 

0.78 



Sampling Frequencies 
(fs= 10 6 l252.xo Hz) 



640 
1280 
1920 
2560 
3200 
3840 
4480 
5120 



15 



20 



25 



30 



To obtain a uniform spaced sampling rate the exact delays given in Table I are approximated as an integer 
number of sampling intervals. The approximated delays, di, are chosen such that the quotient of the 
approximated delay di when divided by the total number of sample-sets n=7 has a different integer 
remainder for each value of i. Table III shows how di is chosen so that its remainder is different for each value 
of i. 

35 Table III 



30 



35 



/ 


Time Delay 


Time Delay 


di 


Remainder of 




(sec) 


1.24 \lS. 




dill 


40 0 


0 


0 


0 


0 


1 


6.21 


5.008 


5 


5 


2 


24.67 


19.89 


20 


6 


3 


54.80 


44.19 


44 


2 


4 


95.69 


77.17 


78 


1 


45 5 


146.10 


117.82 
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6 


204.50 


164.91 


164 


3 
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It should be understood from Table 111 that because di is a function of the time duration of the sampling 
interval, i may be determined for any consecutive sampling interval from the remainder di/7 which is equal 
to the remainder of z/7. Moreover the samples taken during each consecutive sampling interval correspond 
to beam m where: 



50 



55 



60 



m = Remainder 



m=Remainder 



z-di +M 



7^ — M — _ 

z-di +36 

_Z 

_ 3E 1 



+ 1 



+ 1 



where i is known for any di per Table III and di is known for any z (because Remainder [di/7] = Remainder 
[z/7]). 

Referring now to Figure 3 the timing diagram for the sampling schedule may be derived by calculating for 
each sampling interval the sample-set to be sampled and the beam to which it corresponds. For example, 
65 from Figure 3, during the time duration of the 25th sampling interval stave elements S1 9 and S29 ,ir e 



55 



60 



65 
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sampled for beam 24. This information may be calculated from the above equations because for z-25: 

Remainder [di/7] = Remainder [z/7] 
= Remainder [25/7] 

5 

From Table III for a remainder of 4, di= 1 1 6 and i=5. From equation 1 : 
25-116 



m = Remainder 
10 



+36 +1 



36 



= Remainder [23/36 + 1 ] 
= 23+1 
15 =24 



10 



15 



Hence, stave elements m±i are sampled which corresponds to S241655 = S19 and S29. 

Referring to Figure 3 a repetitive sampling sequence can be seen. At time intervals 0, 7, 14, etc. the first 
sample-set for every beam is sampled. The sequence is dependent upon the remainders of di/7. The delay 
20 associated with the sample to be taken corresponds to the remainder of di/7. The first sample in the 20 
sequence is for di/7 having a remainder of 0 corresponding to the sample not delayed in time (i=0). The 
second sample is for remainder di/7=1 corresponding to i=4, the third for remainder di/7=2 corresponding 
to i=3 and so on until the last sample in the sequence is for remainder di/7=6 corresponding to i=2. The 
delays that correspond to these samples are a function of the rounding of the time delays to be a multiple of 
25 the time duration for the sample interval. It should be understood that this sequence repeats itself after every 25 
7 sampling intervals because the number of sample-sets n has been chosen equal to 7. 

Figure 4 is a schematic diagram for the digital beamformer 38 of Figure 1 . The timing and control circuitry 
is shown within line 70 comprising clock 72, counter/divider 74 and a programmable read only memory 
(PROM) 76. The timing diagram of Figure 3 is stored in PROM 76. PROM 76 accesses multiplexers 40 of the 
30 beamformer to sample input lines 82 in accordance with the interleaved sampling method of Figure 3. Lines 30 
82 separate to feed both multiplexers 40. There are 36 input lines 82 (one for each stave element). The two 
multiplexers are used so that symmetrically disposed staves may be simultaneously sampled. 

The sonar signals entering the beamformer 38 via input lines 82 have had their background noise level 
normalized by preamplifiers 32 and have had their bandwidth reduced to 800 Hx by anti-aliasing filter 34 (see 
35 Figure 1 ). Using the sampling frequency of 3200 Hz the combined sampling rate of input lines 82 by 35 
multiplexers 40 will be (252 samples x 3200 Hz =) 806,400 samples per second. 

Multiplexers 40 interrogate sonar signals on input lines 82. The two symmetrically disposed input lines are 
interrogated by multiplexers 40 and summed by analogue summer 80 prior to A/D convertion. PROM 76 will 
access only one multiplexer when only one input line is to be sampled such as for the first sample in the 
40 formation of any one beam. In this instance, the output of this one multiplexer is summed with a nil output 40 
from the other multiplexer resulting in a sum which is indicative of the sample itself. The summation is done 
in analogue because it is quicker and more economical summing the symmetrically disposed staves prior to 
A/D conversion. The analogue sum leaves summer 80 and is converted to digital format by A/D converter 84. 
A pulse generated by PROM 76 on line 86 initiates the A/D conversion. The output of the A/D conversion is 
45 shown to be an 8-bit word which is multiplied by weighing factors supplied by weight circuitry 86 to 45 
multiplier 88. The timing and control issues pulses from PROM 76 along lines 90and 92 to initiate the 
multiplication. The weighing factors are chosen with respect to which samples of the beam are being taken. 
For the first sample the weighing factor will be unity. The other weighing factors may be a fraction less than 
unity so as to aid in suppressing side loops. The output of the multiplier is fed into latch 94 which transfers 
50 the 8-bit word representing the sampled signal to digital summer 96. 50 = 

Digital summer 96 adds all weighted samples to all other samples of their respective sample-set which 
when summed forms a beam output signal. Summer 96 utilizes Random Access Memory (RAM) 98 which 
stores the partial sums of each beam output signal to be formed. Pulses sent from PROM 76 to along line 100 
initiates the summation of each new sample entering input A with the partial sum of its sample-set entering 
55 at input B. This new partial sum A+B is then stored in RAM 98 at an address indicated by bus 102 from PROM 55 
76. Bus 102 also addresses RAM 98 to bring the partial sum of th sample-set to be summed in summer 96 
with each of its new samples. 

RAM 98 stores a partial sum for each beam to be formed. In this case 36 beams are being formed so RAM 
98 muststor 36 partial sums. RAM 98 outputs on bus 104 digital words representing each full beam output 
60 signal. Each full beam output signal is fed from the beamformer only after all 7 samples of the sample-set for 60 
the beam output signal have been summed. At this time the RAM 98 clears that memory location which is 
used to keep the partial sum so that a new output signal can be formed in that direction. 

The beamformer 38 outputs a beam every (7 samples x 1 .24 fis (the time duration of the sampling 
interval)] = 8.68 fis. Hence, the beamformer which samples at a rat of 806,400 samples per second will 
65 output approximately 1 1 5,200 full beam output signals every second. 65 
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It should be apparent that the beamformer requires the fastest operating ^™*£ oceS * na SeCt '° n 25 
(Figure 1) need not work at such a fast rate which makes for more economical hardwa re. 

Besidesthefu.lbeamoutputsigna^^ 
into adaptive time varied gam control (A.T.V.G.) 36 (Figure 1 h l ne a ,^ o c d ve 5 

5 reverberation decay as a function of range, variations in the ^■"J^^^^™^, de cay curve 

Reverberation is the energy from the transmrlted pn « b»no « looted ^."^^"X reflection and 10 
,0 reflection i, caused by discrete particles «" d .^ 1 *7'*"^ m ( ^;S."S^ dynamic range of 

SeSc^^^ 

b tr;%Tc p „n^^^ 

WlZ ySEa aBBfeat. 

^^^^^ 

=^ 30 
During the period between gam ^J^^^^^^^^on computer 1 18 averages 

7th^^^ 

35 correc^n is required" If AG represents the gain correction then the gain corrects requ.red ». 



25 



30 



15 



20 



....G = -20 log Av/R where 

Av = average of the beamformer full beam output signals 40 
40 R = the desired full beam output. „ m „„t 0 riifir-nmrjutes the value of the gain 

For each period between gain changes gain 'correction T .V.G. control 114 

transmitted and received, the new gam values are appl .ed to ' ^Je preampl rtersl n P 

on subsequent ping-cycles so that the gam-time profile « "^""■"^jS of 30 second Because the 

characteristics of the medium. The length of a ping-cycle may b , e J " t ^? h ° 2%° m d1 if rs as the time of the 

range of the target corresponds to time, different gains are applied to the preamplifiers tne m 50 

^^seXh^n^ 

normalization will occur, which are removed in a beam oriented A.G.C. (adapt.ve gam control) g 

(Fiqure 1 ) following the beamformer 38. r *« r thp nain rorrection 

The exact computation of AG is a complicated and time consuming proces .for he g »^J™^ 5 5 
55 computer because of the logarithmic function. The gain ^^ n ^^'^*^^^ 

sample algorithm. For any number of dB's required ^---^^ is to be 

^dtS^ 
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Av/R 

Av 

A-R 



2-5 
2.5 R 
1.5R 



5 for R chosen to be unity. 
Av-1 = 1.5 

The converse of this would apply and if Av-1 -1 .5 a gain correction of -8 dB would be required. Table IV 
10 shows a binary code which has been generated for gain decreases. 



Table IV 



AG (dB) 
15 -20/ogAvlR 



Ratio AvIR 



20 



25 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
30 13 
14 
15 
16 
17 
35 18 
19 
20 
21 
22 
40 23 
24 



Exact 

1.122 
1.259 
1.412 
1.585 
1.778 
1.995 
2.239 
2.512 
2.818 
3.162 
3.548 
3.981 
4.467 
5.012 
5.623 
6.310 
7.079 
7.943 
8.913 
10.000 
11.22 
12.59 
14.13 
15.85 



Fraction 

1 1/8 
1 1/4 
1 3/8 
1 1/2 

1 3/4 
2 

2 1/4 
2 1/2 

2 3/4 

3 1/8 
3 1/2 
4 

41/2 
5 

5 5/8 
61/4 

7 1/8 
8 

8 7/8 
10 

11 1/4 

12 5/8 
141/8 
15 7/8 



Difference Av-R 
Fraction Binary Code 



1/8 
1/4 
3/8 
1/2 
3/4 

1 

1 1/4 
1 1/2 

1 3/4 

2 1/8 

2 1/2 
3 

3 V2 
4 

4 5/8 

5 1/4 

6 1/8 , 
7 

7 7/8 
9 

10 1/4 

11 5/8 

13 1/8 

14 7/8 



0000.001 
0000.010 
0000.011 
0000.100 
0000.110 
0001.000 
0001.010 
0001.100 
0001.110 
0010.001 
0010.100 
0011.000 
0011.100 
0100.000 
0100.101 
0101.010 
0110.001 
0111.000 
0111.111 
1001.000 
1010.010 
1011.101 
1101.001 
1110.111 



10 



15 



20 



25 



30 



35 



40 



It should become apparent from Table IV that for any number of dB's the gain correction may be obtained 
from a read only memory which may be addressed by the binary code obtained from the difference of Av-R. 

45 Figure 7 is a block diagram showing the A.T.V.G. control 36 in greater detail. Beam output signals from the 
beamformer are fed back via bus 1 20 to averager 1 26. Averager 1 26 accumulates a large sample of beam 
outputs and averages them once every period between gain changes. It should be understood that during 
one scan 36 beam outputs are formed. Because it has been found that the minimum period for gam change 
is about 20 mS 32 beam output signals may be added for 32 scans. The four beams aft of the ship are 

50 excluded due to high propellor noise. Hence for each period between gain changes averager 1 26 sums 32 
beam output signals 32 times and divides this summation by 1024. Division in the case of a digital s.gnal only 
requires shifting. The average is fed into comparator 128 along bus 130. Comparator 128 substracts desired 
value R on bus 132 from the average value. Comparator 1 28 outputs the sign value and difference value on 
respective buses 134 and 136 to line selector 138. The line selector 138 uses the sign bit to determine if the 

55 gain should be increased or decreased. If the gain should be increased line select 138 accesses read only 
memory (ROM) 140 by increase bus 1 24 causing ROM 140 to output the gan correction increase at the 
address location determined by the difference value. If the gain should be decreased line select 1 38 accesses 
ROM 140 by decrease bus 144 to output the gain correction decreases from ROM 140 at the address location 
determined by the difference value. 

60 ROM 140 feeds the gain correction along bus 124 and into adder 148. Adder 148 take the gain correction 
value and adds it to the previous gain stored in Random Access Memory (RAM) 1 50 which enters adder 148 
via bus 152 The sum represents the new gain and is outputted on bus 154 on the line selector 1 56. The line 
selector inputs the new gain into the RAM 1 50 via bus 1 58. During the initial ping-cycie the assumed 
reverberation curve information and initial gain is deposited in RAM 1 50 and line selector 1 56 by way of bus 

65 160 RAM 150 outputs th gain to the gain decoder 16 along bus 152. The gain decoder 161 codes the new 



45 



50 



55 



60 



65 
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25 



gain and outputs the coded gain on plurality of lines 162 of control channel 1 12 which adjusts the gain of the 
preamplifiers for each period between gain changes. Hence the T.V.G. control 1 14 computes the new gain 
while simultaneously using the previous gain to control the preamplifiers. 

Now referring again to Figure 1 the first stage of processing section 25 is the automatic gain control (AGO 
5 stage 46. The AGC controller 46 increases or decreases the gain of each full beam output signal independent 5 
of the other beam output signals after a comparison of the actual beam output signal's gain with the average 
gain of all the beam signals of that beam. The AGC controller 46 eliminates any beam-to-beam variations in 
normalization which may occur from the A.T.V.G. control 36. The beam output signal from AGC 46 is fed into 
the own doppler nullifier {O.D.N. ) and fine bearing computer 54, 

10 Although the function of the O.D.N. 48 (Figure 1 ) has been previously described to reduce the full beam 10 
output signal to half by taking away the own ship's doppler effect on the signal, some of the principles of the 
O.D.N, should now be discussed. 

Sonar echo signals, the signals to be detected, are generally received at a frequency which differs from the 
transmitted frequency. The difference is known as the doppler effect which is due for the most part to the 

1 5 relative motion of the target with respect to the ship. 15 
In order to extract the target doppler from the total doppler it is necessary to remove the doppler effect of 
the own ship's motion. The O.D.N, involves the mixing of the received beam formation's frequency with a 
local frequency which corresponds to the ship's doppler followed by filtering out the sum frequencies 
obtained when the frequencies are mixed. 

20 In order to reconstruct the original target informaton from the sonar signals entering the receiver the 20 
sampling rate of the receiver must be at least twice the highest signal frequency present. However, because 
the receiver of the present application incorporates anti-aliasing filters at the input of the receiver, the 
sampling frequency is reduced to at least twice the band of frequencies passed through the anti-aliasing 
filters. It should be understood that the desired frequency band may be "folded" into the baseband. The so 

25 called folding frequencies occur at multiples of half the sampling frequency. Signals occurring at multiples 
outside the desired band will also be folded into the baseband. Such frequencies are called "aliases". The 
purpose of the anti-aliasing filters is to restrict the frequency range to the desired band and thus suppress 
signals in any other band. The mirror image of the baseband extending from 0 Hz into the negative 
frequency region cannot be eliminated. 

30 The doppler effect of the beam output signal may vary from a minimum of f t -f td to a maximum of f t +f t d* 30 
where f t represents the transmitted frequency and f td represents the total doppler shift. 

If the sampling frequency is chosen for the minimum received frequency to be a multiple of the sampling 
frequency, as is the case in the receiver of this applicaton, the received frequency of the beam output signal 
will always fall in the band of (a)(fs) to (where a is an integer and fs the sampling frequency). The 

35 sampling frequency will be folded into the baseband from 0 to fs/2. The frequency of the received beam 35 
output signals is then: 

f e = fdo + f t d, where 

40 f do is the frequency of zero total doppler and f td is the total doppler shift which may have positive as well as 40 
negative values. 

As indicated earlier the maximum doppler shift of the frequency about the transmitted frequency was 
chosen to be 400 Hz. The value of maximum doppler was chosen by assuming that the maximum ship speed 
and target speed were each 40 knots. Given a converstion factor of 5 Hz/knt, the maximum speed of 80 knts 
45 becomes a doppler shift of 400 Hz (= 80 knt • 5 Hz/knt). 45 
If this low frequency received beam output signal, f e , is mixed with a frequency proportional to the own 
ship's doppler (a known variable), the resulting difference component will no longer contain the own ship's 
doppler. This resultant frequency is proportional to the target doppler. 

During the mixing process within the O.D.N, device the following should be considered: 
50 1 . In addition to the difference f requency the output of the mixer will contain the sum frequency, which is 50 
to be filtered out. 

2. Since the doppler is variable, the sum and difference frequencies will each require a frequency band. 
These frequency bands must be separated sufficiently so that the above filtering in the above paragraph may 
be done. 

55 3. The effective noise bandwidth of the anti-aliasing filter may be assumed to be not more than 1 .5 times 55 
the passband. This means that the noise bandwidth of the difference signal will be considerably larger than 
the band required for the actual signal frequencies. In order to obtain a maximum S/N enhancement the 
center of the difference frequency band must be shifted to the right so that the portion of the noise band in 
the negative frequency region, when folded over, falls outside the desired difference band. This is clearly 

60 shown in Figures 8A and 8B. Figure 8A shows the effect fold over of noise on the desired difference 60 
frequency band 170. The filter response is designated as 172 having a noise band 174. As shown in Figure 8A 
the noise band in the negative frequency region 1 76 folds over into the desired band as indicated by shaded 
area 1 78. This folding affect effectively increases the noise in the desired difference frequency band 1 70. 
Referring now to Figure 8B shifting the center frequency of the difference band 1 70 to the right by desired 

65 amount 180 will result in the folded over noise from the negative frequency region, as indicated by shaded 65 
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area 178, falling outside of the desired difference frequency band 170. 

4. The requir ments of paragraphs 2 and 3 above may be satisfied by introducing an "offset" in the 
relation between the beam output signal frequency and sampling frequency, and also in the local oscaillator 
frequency. 

5 The offset in the beam output signal frequency entering the O.D.N, device may be defined as: 5 

A = (transmitted frequency) - (nearest lower multiple of sampling frequency with respect to the 

transmitted frequency) 
and the local oscillator frequency is defined by: 
10 f m = B 4- own ship's doppler, where B is the offset in the oscillator frequency. ' 10 

It should be understood that the optimum solution is obtained when: 

A + B = f s/2 

15 and the center of the difference band (A-B) is located as shown in Figure 8B. 15 
Referring now to Figure 9 there is shown an O.D.N, device which takes into consideration the frequency 
techniques previously mentioned in the mixing operation. 

A read only memory (ROM) 182 is shown addressed by three buses 184, 186 and 188. ROM 182 outputs a 
value K on bus 1 90 which corresponds to the value of: 
20 20 

K Ji_ COS p, where 

c ^ 

ft — transmitted frequency, 

COS p = beam angle relative to the bow of the ship, and 
25 c = the velocity of sound in water (V.O. S.I. W.). 25 

The ROM 1 82 is used to store predetermined values of K for different values of f t , c and COS p, thereby 
eliminating lengthly computation times. The valve of K read from ROM 182 corresponds to the information 
on the three address buses. Address bus 184 is indicative of the transmitting frequency which is a known 

30 variable. 30 
Address bus 186 is indicative of the V.O.S.I.W". for which certain values have been calculated. Lastly, 
address bus 188 is indicative of the beam angle with respect to the bow of the ship, p. Because the direction 
of the beam being processed is known, a signal indicative of the beam may address ROM 182. In summary, 
ROM 182 outputs a value K which represents a predetermined computation for selected known input values 

35 of the transmitting frequency, beam angle and V.O.S.I.W. 35 
The value K is multiplied with the own ship's speed, a known variable U inputted on bus 192 to multiplier 
194. The result of the multiplication is outputted on bus 196 and corresponds to the doppler effect of the ship, 
fos/ where: 

40 f os =Ux^COSp 40 



The own ship's doppler, f os , is then added in adder 200 with offset B introduced on bus 198. The adder 200 
may perform the addition in 2's compliment arithmetic so that the output on bus 202 is f m =B±f os . The plus or 
45 minus sign are indicative of the ship moving forward or backward, respectively. At this point (bus 202) the 45 
O.D.N, device has determined the proper value for the own ship's doppler effect and offset this value by B to 
give the offset in the local oscillator frequency, f m . 

The O.D.N, has now to produce for each beam output signal a string of samples which represent the 
function: 

50 50 r 

COS(2 jt * f m - b • At), 

where t = l/f s . . . (the sampling frequency, the interval between samples of the same beam), and 
b = is a succession of integer values 1 , 2, 3 etc. 
55 This cosine value is to be multiplied with the beam information to obtain the desired difference band. 55 
It should be understood that 2 .t • At is a constant which leaves the quantity f m • b proportional to the 
required angle COS(2 x - f m * b - At). It should be understood that th required angle may be determined by 
the accumulation of values for f m above for each successive beam output signal by using a Random Access 
Memory. 

60 The accumulation is done in angle accumulator or adder 204. Each successive valu of f m is added to the 60 
sum of the prior values stored in RAM 206. The proper correlation between f m and is corresponding 
accumulation is governed by beam address bus input 208 which addresses RAM 206. 

In the system of the pr sent application, thirty-six beam output signals are formed. Hence, adder 204 and 
RAM 206 keep a running value of the accumulated angle, a, for each beam output signal. The value of the 

65 accumulated angle is given as: 65 
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v = b • f which is outputted on bus 210 into a cosine ROM 212. For each value of the >^ ul ^angla 
th re is JcoTrespondingpredetermined va.ue of the function COS.2 ^.^^^^ * 
cosine ROM 212. For the value of the accumulated angle addressing COS ROM 212 on bus <nu me 

5 c =?or u ,^ 

S ' tSSouVpuS", ££££L* «he O.O.N. along bu, ^ The*, n, outpu-ig _n. I information 
,n Lut Ld on buV?20 It should be recalled that offset A corresponds to a frequency value equivalent to the 



15 

224. 



The output of multiplier 224 on bus 226 given as: 
20 COS(A+f, g +fos) • COS(B+fos) 

(For simplicity only positive signs are used.) 

= 1/2 COS(A+B+f tg +2f os ) + 1/2 C0S((A-B) + (f t g+f os ) - fos) 
= 1/2 COS<A+B+f, g +2f os ) + 1/2 COS((A-B) + f« g ). 

the target doppler. in this situation COS(f„l This na it Beam o p a * a|1 ^ 

30dopp,erin,orrna,ion,e^^^ 

a ~ 8 :nX^R g :r?"rhe\i 9 n^ processor50 isshown having a correlator 51 and, sp.c.rum 

'^e wLple behind the correlator is . ha, the best correlation between the •°'9«'" , °<™'i°"; '* d set 

bearing accuracy to about 1 to 1 1/2 . oinnat orocessor tarqet speed is fed 

52 Processor 52 outputs in digital format information to the d.splay correspondmg to. 
1. Target range, 
- 50 2. Target bearing with respect to true north, 

3. Doppler (target speed), 

4. Fine bearing, and 

' fn sum^ry'complete sonar receiver has been described which subsequent digital beamforming 

55 analyses the received sonar signals digitally. 

CLAIMS 

1 An interleaved method of sampling a multi-element symmetrical array in order to form a beam related 
60 to one or more Xectiona^lements wherein the total number of samples to be taken ,s g.ven as. 

Z = Mn 

where M is the total number of beams to be formed and n is the number of sample-sets desired to form any 
65 beam m, said method comprising the steps of: 
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A) sampling each i element symmetrically disposed to the left and right of the one or more directional 
elements for beam m during each consecutive sampling interval z where i is defined as o^n - 1 z is 
defined as o*=z*sZ, and m is given by: 



m=Remainder 



z-di 
n 



M 



+ 1 



where di is the time delay of the i th sample to the right and left of beam m with respect to its first sample 
1 0 expressed as an mteger number of sampling intervals z and di is chosen such that Remainder [di/n] has a 
different integer value for.each value of i, and i is determined from the remainder [z/nj- 

, sumi ™9 together each sample taken from i elements to the right and left of beam m to form a 
sample-set; 

C) adding each sample-set taken to its corresponding partial sum of sample-sets and storing the partial 
lb sum; and 

D) forming the beam from the partial sum of sample-sets for each beam when the beam's respective 
sample-sets have been summed. 

2. A method of beamforming a plurality of digital beam output signals in sonar receiver having a 
multi-element symmetrical transducer array wherein each beam output signal has a direction related to at 
20 least one respective directional element in the transducer and the total number of samples to be taken is 
given as: 

Z=Mn 

25 where M is the total number of beams to be formed and n is the number of sample-sets desired to form any 
oeam m, said method comprising the steps of: 

A) sampling each i element symmetrically disposed to the left and right of the at least one directional 
element for beam m during each consecutive sampling interval, z, where i is defined as O^i^n - 1 z is 
defined as 0*£z=SZ, and m is given by • 

30 

'"z-di 

-U. IWI 



m=Remainder 



+ M 



M 



35 where d. is the time delay of the i th sample to the right and left of beam m with respect to its first sample 
expressed as an integer number of sampling intervals z and di is chosen such the Remainder [di/n] has a 
d.fferent integer value for each value of i, and i is determined from the Remainder [di/n] which remainder is 
equal to Remainder [z/n]; 

B) summing together each sample taken from i elements to the right and left of beam m to form a 
40 sample-set; 

C) converting the sample-set from analogue to digital format; 

D) multiplying the digital sample-set with weighting factors; 

E) adding each digital sample-set to its corresponding partial sum of digital sample-sets and storing this 
partial sum; and . 9 

F) forming the beam from the partial sum from each beam when all n sample-sets for beam m have been 
sum iii eo. 

3. The method as claimed in claim 2 wherein the step of analogue to digital conversion is done prior to 
the step of summing together the samples taken from i elements to the right or left of beam m to form the 

sample-sets, 

50 4. A digital beamformer for a sonar receiver having a multi-element symmetrical transducer array 

InT^f °, ? ° U ? Ut Si9nal h3S 3 direction re,ated to ^ least one directional element in the transducer 
and the total number of samples taken by the beamformer is given as: 



45 



Z = Mn 



55 



where M is the total number of beams to be formed and n is the number of sample-sets desired to form any 
Deam m, said beamformer comprising : 

control means for sampling each i element symmetrically disposed to the left and right of the at least one 
direction element for beam m during each consecutive sampling interval, z, where i is defined as 0^n - 1 z 
60 is defined as 0^z*=Z, and mis determined by: 



m=Remainder 



z-di 



M 



+ 1 
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S ^^U-a m p l es. ate „ ( ,o m ea eh i e . em e n „o t h eriaM8 nd l e«o,be 8m n 1 ,o f o TO a 5 

Tn^D converter for convening the sample-set from analogue, o digita, format; 
armltipli«)brWBlohtlng*edWtf 

,0 ^^^^^^^^^^^ 
10 SS£5K, Mr "^^T^ESS, ft- > -*■> »«- output signal o, a 
5. A method of extracting own ship f^°J^.^ juT me thod comprising the steps of: 

''^ge^g^^^ 

digital signal by frequency amount B wn.cn is a run 

'TmSngthegenerateddig^— 2 0 
20 difference frequency component and a surr frequen cy o™ di „ erencs „ equencv ban d has a center 

^cyS^ 

A = ft-M fs .. 

anriM is th e nearest lower multiple ofthebeamformer samplmg 

30 wherein f t is the transmitted frequency, and M fs .s the 

frequency to the transmitted ^ ue ^. ■ offset » B » is cho sen so that 
7. A method according to claim 6, wherein onset 

35 B=£-A 

where f, is the beamf ormer sampling f r «J^ n «£ r informat ; 0 n from a digital beam output signal of a 
Mb eUrm=^» 

i«se^h°ege^ 

45 T^rS^ fd^ 

obtain a signal having difference ,h. difference freouency band is offset ,n 

CJezr~^^^^ 

component representative of th. targe t doppler ,n orm««on ; ^ b8twee „ tr a nsmi ,,ed 

■ (re yrar.^ 

55 frequency. . for gan .,a,ing a digital signal representative of theown 

"7,. Adaptive time varied gai r ^'"'.■'^SS.lv proportional to momentary average 

60 reverberalioninthemediumofoperahon ,, h J„ gaft0 re4^n; 

storage means having a ^.lit, »fslo.W 2» » approximations of an mverse standard 

eS^e'Sorstoringeachorsaid^ 
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storage space where said time-duration's repective gain step approximation is initall v stored ■ 

<JL ,I he apparatus of claim 14 wherein the determining means further includes a memory which has 
20 5^S!rS, n „! re,totern,,ned Va ' UeS 93in COrreCti ° n Which are aCC6Ssed "V the coded v W at h o h f the 



